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Bound states around an impurity are investigated for a two dimensional electron system in a strong 
magnetic field. Long-range Coulomb potential and related potentials are considered. Schrodinger 
equation is solved numerically to obtain the bound states. The energy and wave function of these 
bound states are indirectly observed by the scanning tunneling spectroscopy as local density of 
states (LDOS). Theoretically obtained LDOS is compared with experiment. Reasonable agreement 
is obtained. 



I. INTRODUCTION 

In recent experiments, Matsui et aL 1 and Niimi et 
al£ succeeded to observe local density of states at a 
surface of graphite by scanning tunneling spectroscopy 
(STS). In a strong magnetic field at low temperature, 
they observed a central peak and ring-like structure 
around impurities. Graphite is a three dimensional ma- 
terial; it is a semimetal having both electron band and 
hole band. However, HOPG (highly oriented pyrolytic 
graphite) that is used in the experiments looks like two- 
dimensional when observed by the transport and STS; 
the two dimensional Landau levels are observed, and 
quantum Hall like plateau is observed^ The origin of the 
two-dimensionality is attributed to abundance of stack- 
ing faults in this type of graphite. In this paper we exam- 
ine to what extent the observed structure can be under- 
stood by an ideal two-dimensional model, and clarify the 
meaning of the structure observed in these experiments. 

Although graphite has both electron band and hole 
band, here we consider simple 2-dimensional electron sys- 
tem with parabolic dispersion. We place a positively 
charged impurity at the origin of the two-dimensional 
plane, or some distance d separated from the origin in 
the direction perpendicular to the 2-d plane, taking into 
account the possibility that the impurity is situated in the 
bulk of graphite. In this model, electronic states around 
the origin in a strong magnetic field are investigated. 

The angular momentum around the origin is a good 
quantum number, and the bound states of definite an- 
gular momentum are formed below each Landau levels. 
Which value of the angular momentum gives the lowest 
energy state for each Landau level depends on the dis- 
tance d. When d is less than the order of the magnetic 
length, the states with angular momentum zero are the 
lowest energy states for each Landau level. If d is much 
smaller than the magnetic length, the wave function of 
these states has a sharp peak at the origin, and the en- 
ergy is separated from other states. In such a case the 
local density of states (LDOS) at the energy between the 
Landau levels shows a sharp central peak surrounded by 
ring-like structure. This shape of LDOS resembles ex- 



perimental results. On the other hand, when d becomes 
large, such structure as observed experimentally is lost. 
Thus, we claim that the observation of the LDOS by STS 
can give information on the possible structure of the im- 
purity potential. 

This investigation has been begun as collaboration 
with Fukuyama's group for the explanation of their ex- 
perimental data. Part of the present results arc included 
in Niimi et al's letter^ This paper gives general treat- 
ment of the problem and details of the calculation. 

Recently, truly two-dimensional graphite is realized, 
which is known as graphene4^& In this case the energy 
spectrum is not parabolic, but linear, and carriers there 
are sometimes called massless Dirac fermions. How the 
LDOS around an impurity in this material shows up is 
also an interesting problem. Theoretical investigation for 
this case is given in a forth coming paper. 

This paper is organized as follows. In §2 potential 
used in this calculation and the method are explained. 
In §3 we notice that the wave function and energy of 
positive angular momentum state is known from those 
of the negative angular momentum state in the presence 
of isotropic potential of arbitrary shape. The concrete 
energy and wave function in the presence of impurity po- 
tential are obtained in §4 and 5. We calculate local den- 
sity of states in §6 using the results of energy and wave 
functions. Finally in §7 comparison with experiment is 
done. 

II. MODEL AND METHOD 

We consider a model of two-dimensional electron sys- 
tem in a magnetic field with a positively charged impurity 
at the origin. The Coulomb interaction between electrons 
is neglected, so we consider single electron states in this 
paper. Being a single electron problem, the spin freedom 
is also neglected. Then the Hamiltonian is 

H = H Q + V(f) , (1) 

where 

Ho = 7 ±-(P-eAf, (2) 



is the unperturbed Hamiltonian, and V(f) is the impurity 
potential. Here, m c is electron's band mass, and e < is 
its charge. The momentum p = (p x ,p y ), the coordinate 
r = (x,y), and the vector potential A = (A x ,A y ) are 
two-dimensional vectors. A uniform magnetic field B is 
applied perpendicular to the 2-d plane, so the vector po- 
tential is given as A — (—By/2, Bx/2) in the symmetric 
gauge. 

Considering a possibility that the impurity is situated 
at a distance d from the 2-d plane, we investigate a sys- 
tem where the impurity potential is given as 

V(r) = -- (3) 
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where e is the dielectric constant. This potential includes 
two limiting cases; one is pure Coulomb case, d = 0, 
and the other is a case where d is large enough that the 
potential around the origin can be approximated by a 
harmonic potential: 



V{r) 



1 



m c Wnr 2 + const. 



In terms of these variables, the Hamiltonian Hq and the 
angular momentum L z are written in quadratic forms, 



(10) 



and 



Then by introducing ladder operators for the Landau 
level, 



h 



V2h 
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(tt x -m y ) , 



(4) and those for the angular momentum, 
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where u 2 . — e 2 /4ired 3 m c , and ujq gives frequency of the 
harmonic oscillation of an electron around the origin in 
the absence of the magnetic field. 

The eigenstates of this Hamiltonian in the limit of 
large d, the harmonic potential case, has been analyti- 
cally obtained and well known^ but those with general 
values of d are not. Here, we solve this problem numer- 
ically. Namely, since we know the eigenstates of the un- 
perturbed Hamiltonian, the Landau states, we calculate 
the matrix elements of the potential, and diagonalize the 
Hamiltonian with the basis of Landau states. 

The procedure is straight forward. However, before 
giving the results, we summarize the eigenstates of the 
unperturbed Hamiltonian Hq, and also the eigenstates 
with the harmonic potential. One way to diagonalize 
Hq is to use a set of canonical variables, the dynami- 
cal momentum 7? and the center coordinate R — (X, Y) 
instead of canonical momentum p = (p x ,p y ) and coordi- 
nate r = (x,y). They are defined as 



7? = p — eA 



(5) 



V2h 



(X-iY) 



h = 7% {x + ' lY) 

we rewrite Hq and L z as follows 
Hq = huj c ( a) a + 



and 



L z = h (a f a - 6+6) . 



(14) 



(15) 



(16) 



(17) 



Here a> c = \e\B/m e is the cyclotron frequency. Simulta- 
neous eigenstates of a^a and b%, \n,m), which satisfies 



a} a\n, m) — n\n, m) 



(18) 



and 



X = X - — TTy , 

Y = y+ l fn x , 



(6) 
(7) 



where Ib = \/Ti/\e\B is the magnetic length. These vari- 
ables satisfy commutation relations, 



Fa, TTy J = -1 



I 2 ' 
'B 



(8) 



6 T 6|n, m) — m\n, m) , 



(19) 



is also the eigenstates of Ho and L z , and eigenvalues are 
given as E n = (n + l/2)?kj c and l z fi = (n — m)Ti, respec- 
tively with n and m being non-negative integers. 

The eigenstates in the presence of the harmonic poten- 
tial V(r) = (l/2)m c u>Qr 2 are obtained by our using the 
ladder operators a, w , and b at larger magnetic field, 
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In terms of them, H and L z are expressed as follows 

H = H + ^m c uj 2 r 2 

= (tfa + b r b + l) + t^c (a f a - & f &) , 

(21) 

and 

L z = h (a? a - tfb) . (22) 

Thus, the eigenstate of the ladder operators at B, 
\n,m)g, is the eigenstates of H and L z at B with eigen- 
values, 

H\n,m)g = E niTn \n,m)g 

= j| \A>c + 4w o( n + 171 + !) + T^cfa - m) | |n,m) 5 , 

(23) 

and 

L z \n, m) g = h(n — m)\n, m) g . (24) 
I 



It should be noted that in this case the lowest energy 
state for each Landau level is the state with m = 0, or 
angular momentum l z — n. 

In order to obtain the eigenstates of the Hamiltonian in 
the presence of the Coulomb potential V(r), we express 
the total Hamiltonian as a matrix using the eigenstates of 
H and L z , and numerically diagonalize it. For numerical 
calculation we parameterize the Coulomb potential: 



v{r) = HlJc vFTp ' 



where £ = r/l B , 5 = d/Z B and a = e 2 / ATrehuijB- Namely 
we use w c and Ib as units of energy and length, respec- 
tively. In the following we use a symbol £ for energy 
expressed in the units of fko c , namely, E = £fru) c . 

Even in the presence of V(r), the angular momentum 
remains to be a good quantum number. So the Hamilto- 
nian is block diagonalized for each angular momentum. 
For a given l z = n — m, it is expressed as 



A = £| n , n _ U (n+±) (n,n-, 



+ a V |m,m - Z z )(ni,m - l z \ * — = \n 2 ,n 2 - l z )(n 2 ,n 2 - l z \ . (26) 

V? + 6 



n 1 ,n 2 

For the case of d/ls = 5 = 0, the matrix elements are calculated analytically: 

/ t 1 1 1 , x 1 ( N X \N 2 \ N 1/2 



mini 



( ^" 2) r(i + |/,| + i) 

J=0 



. V- ^UT W t;; , 2JVi - 2 j - 1 W 2N 2 - 2j - 1 » 



where 

JV li2 = ni i2 -^(i z + |i a |), (28) 

and T is the gamma function. We label the eigenstates by 
n and m so that they evolve continuously as a increases. 

III. RELATION BETWEEN POSITIVE AND 
NEGATIVE ANGULAR MOMENTA 

In our system with isotropic impurity potential, the 
angular momentum is conserved. When the direction of 
the magnetic field is in the positive z-direction, which is 



the choice in this paper, a state in the lowest Landau 
level, n = 0, has only zero or negative angular momen- 
tum. In the higher Landau levels (n > 0), most states 
have negative l z , but those with < m < n have posi- 
tive angular momentum. The energy and wave function 
of these states with a positive angular momentum l z are 
known once we know the eigenstate with angular mo- 
mentum — l z in the lower Landau levels. We explain this 
relation in this section. 
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The Hamiltonian in the polar coordinate (r, 9) is 



H 



2to c 
.1. 



i d 



r dr \ 6V 



1^1 

r 2 002 



9 1 



(29) 



We write the wave function as 

ip(r,8) = e il * e R(r). 



(30) 



Then the Schrodinger equation for R(r) with eigenenergy 
E is 



2m r 



R" 



-R' 



tR 



2 2 

-m c uj n r 



V(r) 



R = 



E - -huj c l z 



R. 



(31) 



In the left hand side of this equation only l z appears. 
Therefore, we can see that we obtain the same R(r) for 
both positive and negative angular momenta. For the 
same R, the eigenenergies are different, since the right 
hand side has bare l z . If we write the energy at positive 
l z as E > and that at negative angular momentum I' = 
— l z as E <: these energies are related with the following 
equation: 



E> - ^huj c l z 



— TlLOrX' 

2 



2 



(32) 



Namely, the energy of a state with a positive angular 
momentum l z is higher than the negative angular mo- 
mentum counterpart by I z Tiuj c : E > = E < + l z ftio c . In 
terms of the labelling of the eigenstates with n and m, 
we can say that the wave functions are related as 



(r|n, m}e imti = (r\m,n)e 



UlO 



(33) 




and 



0.6 0.8 

a 



E n 



+ (n — m)Tiuj c 



(34) 



It is easily checked that the eigenstates in harmonic po- 
tential have this property. Considering the fact that the 
classical orbits of these states are quite different, this 
result for the radial part of the wave function is quite 
noteworthy. 



IV. ENERGY LEVELS 

In the absence of the impurity potential, the energy 
spectrum consists of Landau levels, and they are degen- 
erate with respect to m. In the presence of the impurity 
potential, the degeneracy is lifted. The resultant energy 
spectrum for 5 = in the units of ?UjJ c , £, is shown in 
Fig. n as functions of the potential strength a. In this 
figure states with angular momentum l z Ti = (n — m)K = 
—2h to fi are shown. As shown in this figure the lowest 
energy state among a given Landau levels is always the 
l z = state. This is in contrast to the case of harmonic 
potential, in which case the lowest energy states are those 



FIG. 1: Splitting of the Landau levels in the presence of 1/r 
potential. Energy of the states in the units of hu> c , £, is shown 
as functions of the potential strength a. States with angular 
momentum l z h = (n — m)h = — 2% to h are shown. 



with m = 0, or l z = n. Thus, except for the lowest Lan- 
dau level, the sequence of the energy levels is different. 

The reason for the different behavior is understood 
when we consider corresponding classical orbits for each 
n and to. In the absence of the potential, the classical cy- 
clotron orbit is a circle centered at the center coordinate 
R = (X,Y) with radius tq. The energy and the angular 
momentum in the classical mechanics are given as 



and 



E = \m c w 2 c rl , 



L z = ^\e\B(rl-R*) 



(35) 



(36) 
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FIG. 2: Classical cyclotron orbits of an electron which has 
quantum number n — 1 and m — to 3, or l z — —2 to 1. The 
values of l z is written by each orbit. 

Thus, the Landau quantum number n is related to rg, 
and m is related to R 2 . Namely, r 2 — (2n+l)/| and R 2 = 
(2m + 1)1%. The classical orbits for the second lowest 
Landau level (n = 1) with m = to 3 or l z = — 2 to 1 are 
shown in Fig. [3 where the positions of the orbital centers 
are placed at distance R from the origin, the direction 
angle being increasing by n/2. If we consider the effect 
of the impurity potential perturbatively, the average of 
the impurity potential along the orbit gives the energy 
gain from the potential. Since the classical orbit with 
l z = always passes through the origin, where the l/r 
potential diverges, it has the lowest energy. On the other 
hand, the orbit with m = has orbit center nearest to the 
origin, so the average energy becomes lowest for the case 
of r 2 potential. These differences in the classical orbits 
are reflected in the behavior of the quantum-mechanical 
wave function. Namely, wave function with l z = only 
has non-zero value at the origin. The wave functions 
with l z — suffers large deformation in the presence of 
the l/r potential as we will see in the next section. Thus, 
the energy lowers quadratically as a increases. On the 
other hand, those with l z ^ do not have large amplitude 
where the potential has large negative value, and suffers 
less deformation. Thus, the energy lowers almost linearly 
with a. 



We can expect that as the potential minimum at r = 
becomes shallower, the potential approaches the har- 
monic potential, and the energy of the l z = state be- 
comes higher than the m — state. Thus, as the pa- 
rameter d in the potential increases the sequence of the 
energy changes. The behavior of the energy spectrum at 
a = 0.5 as a function of d is shown in Fig-EI As expected, 
the level crossing occurs at around d ~ Zb except for the 
lowest Landau level, where state with m = and that 

with l z = are the same state. In this figure, we notice 
1 1 1 1 

2.0'"" 



e 1.0-- 



. m =2 

m=3 

-1.Q.I ' ' ' ' 

0.0 0.2 0.4 0.6 „ 0.8 1.0 



FIG. 3: Energy spectrum of states around an impurity with 
strength a — 0.5 as functions of d. The energy of the lowest 
states at d — 0, the l z — n — m = states, becomes higher 
as d increases. At around d ~ Ib, the lowest states become 
those with m — 0. 



that only the l z = states suffer large shift of energy as 
d increases. This behavior is understandable if we con- 
sider the classical orbits. Namely, as d increases only the 
potential around the origin changes. Thus, only states 
which passes through a space around the origin suffers 
shift in the energy. 



V. WAVE FUNCTIONS 

In the absence of the impurity potential, the wave func- 
tion is 



(r\n,m) = ^„, m (r, 6) 



(N+\l z \)\ 



exp 



T2 



IM 



J N 



(37) 



where £ = r/?B ; lz = n — m, N = n — (l z + |i z |)/2, with the same l z are mixed, and the wave functions are 
and polar coordinate (r, 9) is used. At non zero a, state 
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FIG. 4: Squared absolute value of the wave function with 
n = m = 0, 27r/| \tpo,o{r, 0)\ 2 at a — 0.5, solid line, and at 
a = 0, dashed line. 
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FIG. 5: Squared absolute value of the wave function with 
n = m = 1, 2irlft\ipi,i(r,8)\ 2 at a = 0.5, solid line, and at 
a — 0, dashed line. 




deformed from this form. We will sec in Figs. H to El how 
these wave functions are deformed in the presence of the 
1/r impurity potential. 

In Fig. 01 the squared absolute values of the wave func- 
tion 1^*0,0 ('"i @)\ 2 at a = and that at a — 0.5 are shown. 
The round Gaussian peak of the wave function at r = 
for a = is raised and deformed to a cusp in the pres- 
ence of the potential. Actually, we can show analytically 
that the wave function with angular momentum zero has 
a cusp at r = by examining the Schrodinger equation 
at small r. By considering the most important terms 
around r = we can show that 

|^„, n (r, 9)\ 2 ~ |t/w(0, 0)| 2 (1 - 4a£) . (38) 

For other examples of the wave functions, we show in 
Fig.Elthc wave function at n — 1 and m = 1, and in Fig.^l 
those at n — and m = 1 or at n = 1 and m = 0. 8 The 
former is another example of the wave function at l z = 0. 
The latter is an example at l z — n — m = ±1 ^ 0. The 
deformation of the latter case is smaller than the l z = 
case, since this state does not have amplitude at r = 
where the potential diverges^ The smaller deformation 
is also reflected in the smaller energy lowering of these 
state at l z ^ 0. 

VI. LOCAL DENSITY OF STATES 

The scanning tunneling spectroscopy (STS) measures 
local density of states as a function of energy E and posi- 
tion f. We here examine how the states around the impu- 
rity potential show up in the STS experiment. We con- 
sider an empty two-dimensional system with the Fermi 
energy Ep — 0. Since, there is no electron in the two- 
dimensional system before electron tunnels into the sys- 
tem, interaction between electrons can be neglected, so 



FIG. 6: Squared absolute value of the wave function with 
n = and m = 1, 27tZ||t/jo,i (r, 8)\ 2 at a = 0.5, solid line, and 
at a = 0, dashed line. The wave function with n = 1 and 
m — has the same form. In this figure the vertical scale is 
enhanced by four times compared to Fig. 0]and|K] 

single electron states considered so far give correct result. 
Assuming a constant level broadening T for states above 
the Fermi energy, we can write the local density of states 
(LDOS) as 

(39) 

where E n , m is the energy of the state with the quantum 
numbers n and m. This LDOS is proportional to the 
differential tunneling conductance, and measured in the 
experiment. 

In the absence of the impurity potential, summation 
over m makes this LDOS independent of the position, 
and energy dependence is given by the superposition of 
the Lorentzian distribution functions. In the presence of 
the 1 / \fr 2 + d 2 potential at the origin, it becomes depen- 
dent on the distance from the origin. We will investigate 
how these dependences reflect the potential form. 

Having in mind comparing the result with existing ex- 
periments, we choose parameters a and r to be a = 0.6 
and r = Q.2Tilu c . As for the distance d we consider two 
cases, 6 = d/l-g, = and 6 = 0.2. First we show LDOS 
integrated in a circular area of radius 1.5 Zb around the 
impurity in Fig. {7\ 

/■1.5/b 

Di(E)= / dr2nrD(E,r) . (40) 
Jo 

The results for S = (solid line) and (5 = 0.2 (dash-dotted 
line) almost coincides except for low energy region, where 
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FIG. 7: Local density of states integrated in a circular area 
of radius 1.5Zb around the origin, where an impurity sits. To 
make the results dimensionless hu c Di(E) is plotted as func- 
tions of £. The dashed line shows the result for a = 0, the 
solid line shows that for a — 0.6, S — 0, and the dash-dotted 
line shows that for a = 0.6, 8 — 0.2. Here the level broadening 
parameter T = 0.2hu> c is used. 



the energy of the n = m = state differs considerably 
between these two cases. The result without the impurity 
potential (a = Q, dashed line) shows superposition of 
the Lorentzian distribution as stated above. The shift of 
the peaks to the lower energy side at finite a indicates 
formation of bound states around the origin. We will 
compare this result with experimental result later. 

To examine the details of the LDOS we plot the value 
of D(E, r) in gray scale in Figs. |H1 and [5| for the case 
of 6 = and 6 = 0.2, respectively. In these figures, 
the horizontal axis shows distance from the origin, the 
vertical axis shows energy in units of fuv Cl and the higher 
value of D(E, r) is shown by brighter pixels. We notice 
that a series of ridges extend horizontally in these figures. 
Each ridge is mainly composed of states of the same n. 
The peak energies of the ridges tend to the energy of the 
Landau levels in the absence of the impurity £ = n+ 1/2 
at £ — > oo. A remarkable difference between the case of 
5 = and S = 0.2 is that the ridges are not continuous 
for the case of S = 0, namely, several peaks are separated 
by saddle points along the ridges. 

These peaks in Fig. [H] reflect the maximum of the 
squared wave functions. Namely, for the case of states 
with Landau index n — 2, the peak at £ = and £ = 2 
reflects the fact that a state with n = 2 and l z = has 
energy of about 2 tko c and the wave function is peaked 
at £ = 0. The next peak at around £ = 1 and £ ~ 2.2 
comes from state with n = 2 and l z = 1 and l z = — 1 
whose wave functions are peaked around £ = 1. For the 
case of n = 1 series, similar correspondence can be seen 
between the peaks in LDOS and wave functions. 



To give more quantitative picture of LDOS, we show 
in Fig.[IJ]]the ^-dependence of D(E, r) at fixed energy E, 
at £ = 2 , 2.1 • • • and 3. Comparing these figures with 
the wave functions in Figs. |S] and HJI we can see that the 
saddle points in Fig. [S] are formed by the deformation 
of the wave function, and by the difference in energies 
of these states. Figure 5 shows that the squared wave 
function of l z — state decreases around 0.5 < £ < 1.0. 
The increase of the squared wave function of l z ^ is 
not large enough to compensate this decrease, since the 
energy difference makes contribution to LDOS of these 
states at the energy of L — state smaller. 

dos a06_3ds.Sfffi 2006/10/17 18:04:43 

On the other hand, when S — 0.2, these peaks of the 
wave functions overlap and saddle points are not so clear. 
This is because the energies of the states are similar, and 
deformation of the wave function around the origin is 
smaller in this case. Thus, even though the integrated 
LDOS loofs °similar~for 0( these two cases, the details of 
LDOS show clear difference. When 5 increases further, 
the ridges become more smooth and they approach hor- 
izontal ridges of system with a — 0. 




FIG. 8: Gray scale plot of the LDOS D(E, r) for a system with 
impurity potential of a = 0.6 and d — 0. The horizontal axis 
shows distance from the origin, the vertical axis shows energy 
in units of hu> c , and the higher value of D(E,f) is shown by 
brighter pixels. Here the level broadening parameter V — 
0.2hcu c . 

A similar figure of LDOS at given energies as Fig. 
for the case of 5 = 0.2 is shown in Fig. ^2 In Figs.llOland 
1111 we can also see how LDOS looks like at intermediate 
energies. Especially, we observe that at the energy of 
the Landau level far from the impurity, £ = 2.5, LDOS 
becomes smaller than surrounding regions. 



VII. DISCUSSION 

We have seen that the LDOS around an impurity po- 
tential has structure reflecting the strength and shape of 
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FIG. 9: Gray scale plot of the LDOS D(E,r) for a system 
with impurity potential of a = 0.6 and d/l = 0.2. The hor- 
izontal axis shows distance from the origin, the vertical axis 
shows energy in units of hu) c , and the higher value of D(E, r) 
is shown by brighter pixels. Here the level broadening param- 
eter T = Q.2huj c . 



the potential. For the case of d = 0, the peaks of LDOS 
have the same form as the wave function. This theoret- 
ical result indicates that under suitable conditions part 
of the wave function is directly observed by the STS ex- 
periment. In recent experiments on graphite, Niimi et 
al. experimentally observed LDOS in a strong magnetic 
fields The observed structure of the LDOS between the 
Landau levels shows quite similar structure as the present 
result for pure Coulomb potential (d = 0). 

Here we make quantitative comparison between the 
theory and their results. First we show experimental 
differential tunnel conductance averaged over an area of 
20x20 nm 2 in Fig. ^| to extract relevant parameters. 
The horizontal axis shows bias voltage of the probe rel- 
ative to the Fermi level. This differential tunnel conduc- 
tance is proportional to the LDOS, so V > part of this 
figure should be compared with Fig.0 Far from the im- 
purity, we see equally separated peaks corresponding to 
Landau levels. The separation of peaks, about 14.6 mev, 
is consistent with the cyclotron energy hu c calculated at 
B — 6 T and with electron effective mass of 0.057 m e >i2iii 
The shape of the peaks looks like Lorentzian, and the 
parameter T is deduced to be about 2 mev. The three- 
dimensionality of the graphite also should give width to 
the peaks, but the shape of the peak in that case is not 
Lorentzian, so we consider that the width comes from the 
life time of the states above the Fermi energy. 

The peaks observed around the impurity are shifted to 
the low energy side. From the size of the shift, we esti- 
mate approximate value for a to be 0.6. If the impurity 
has a unit charge, a = 0.6 at B = 6T is obtained by 
putting the dielectric constant e ~ 15erj. Since graphite 




*■ £= 3.0 



£= 2.0 



1.0 2.0 3.0 4.0 
\ 



FIG. 10: Local density of states at fixed energies are shown 
as functions of £. In this figure a = 0.6 and 5 = 0. Eleven 
values of £ between £ — 2 and 3 with interval 0.1 are chosen. 
Here the level broadening parameter F = 0.2hu c . 



is a semimetal, the dielectric constant is not well defined. 
Sometimes the value of e = 10eo is used in the absence 
of the magnetic field as contribution from valence bands. 
In the present situation, apart from contribution from 
valence bands, existing free carriers in the partially filled 
Landau level may contribute to screening of the impurity 
potential. Virtual transitions between Landau levels also 
contribute to the dielectric constant. The data in Fig. 1121 
show that the Fermi level (V — 0) is a little above the 
hole Landau level peak at V ~ —2 mev. Thus the Lan- 
dau level is almost filled, and low density of free holes 
may be remaining. These holes can screen negatively 
charged impurities, but are ineffective in screening the 
positively charged impurities. Therefore, we can neglect 
the screening by free carriers in the present sample. On 
the other hand, the contribution from the virtual transi- 
tion depends on the wave vector q, and enhancement of 
e at small q is estimated by a simple calculation to heU 



1 + 2 ql B ^ e -^/\ 
as 



(41) 



where ob is the Bohr radius. In the present case the Bohr 
radius and the magnetic length is of the same order, and 
the enhancement factor at q ~ l^ 1 is about 2. Thus the 
value e = 15eo and a to be 0.6 are acceptable. 

Using these parameters, we compare spatial depen- 
dence of the LDOS measured between the first and the 
second Landau levels. In Fig. 1131 experimental differen- 
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tion the actual band structure of the graphite should be 
taken into account. Although this discrepancy remains, 



the qualitative agreement of the central peak shows that 




V[mv] 



FIG. 11: Local density of states at fixed energies are shown 
as functions of £. In this figure a = 0.6 and 8 — 0.2. Eleven 
values of £ between £ — 2 and 3 with interval 0.1 are chosen. 
Here the level broadening parameter F — 0.2hu> c . 

tial conductance and calculated LDOS are compared at 
several magnetic fields. The value of a is chosen to be 
0.6 at B = 6T. Since the a is inversely proportional to 
square root of the magnetic field, it increases as the mag- 
netic field decreases. The values are a = 0.74, 0.85 and 
1.04 at B = 4T, 3T and 2T, respectively. The energy 
at which the theoretical LDOS is calculated is chosen as 
£ = 2.0, 1.9, 1.8 and 1.65 for B = 6 T, 4 T, 3 T and 2T, 
respectively, so that the central peak becomes highest. 
Theoretical LDOS for potential with d/ls = 0.2 are also 
plotted by dashed lines for comparison. In this case, the 
values of £ are 2.1, 2.05,1,95 and 1.85 for B = 6T, 4T, 
3 T and 2 T, respectively. Comparison shows that the 
width of the central peak for 5 = is quite similar to the 
experimental width. The sharpness of the experimental 
peak indicates that d should be quite small; the potential 
at the impurity is well described by 1/r potential. 

On the other hand, the position of the side peak, which 
forms ring-like structure around the central peak, is not 
well reproduced. Experimental peaks are nearer to the 
origin. If the value of a is much increased, the side peak 
comes nearer to the origin. However, it will increase the 
binding energy considerably, and behavior of the aver- 
aged differential conductance cannot be reproduced. Pos- 
sible reason for the nearer side peak is the three dimen- 
sionality of the graphite. We need further investigation 
for the explanation of this discrepancy. In such investiga- 



FIG. 12: Experimental differential tunnel conductance dl /dV 
averaged over an area of 20x20 nm 2 . The data are taken at 
B — 6 T. The horizontal axis shows bias voltage of the probe 
relative to the Fermi level. The open circles show the data far 
from the impurity, and the closed circles show those around 
the impurity. 



ment of STS, and clarifies that the shape of the wave 
function can also be observed. 

In this paper we investigated bound states around an 
impurity for an ordinary two-dimensional electron sys- 
tem, and compared the results with experimental data 
on graphite, which is a three-dimensional system, actu- 
ally. Recently true 2-d graphite is realized and called 
as graphene. In this case the energy spectrum is not 
parabolic but linear. Similar calculation for this case of 
linear spectrum is underway. The result will be published 
in a near future. 
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FIG. 13: Experimentally observed differential tunnel con- 
ductance (closed circles) and theoretical LDOS for a potential 
with 5 = (solid line), and that for a potential with S — 0.2 
(dashed line) as functions of the distance from the impurity 
£ = r/ls- The four values of magnetic fields, 2T, 3T, 4T and 
6T are chosen. 



